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Abstract 

We analyse the CP violation in supersymmetric models with Hermitian Yukawa 
and trilinear couplings. We show that Hermitian Yukawa matrices can be im- 
plemented in supersymmetric models with SU(3) flavor symmetry An impor- 
tant feature of this class of models is that the supersymmetric contributions 
to the EDM of the neutron and mercury atom are suppressed. In this frame- 
work, ek can be saturated by a small non-universality of the soft scalar masses 
through the gluino contribution. We perform a detailed analysis for the super- 
symmetric contributions to e'/e. Although, the gluino contribution is found to 
be negligible due to a severe cancellation between LR and RL mass insertions, 
the chargino contribution can be significant and accommodate the experimen- 
tal results. Additionally, we find that the SUSY contributions to e'/e from 
the effective sdZ vertex and the AI = 3/2 operators are insignificant. We 
point out that the standard model gives the leading contribution to the CP 
asymmetry in B — > ipKs decay, while the dominant chargino contribution to 
this asymmetry is < 0.2. Thus, no constraint is set on the non-universality 
of this model by the recent BaBar and Belle measurements. 



1 Introduction 



CP violation is one of the outstanding problems in high energy physics. Although the 
standard model (SM) is able (till now) to accommodate the experimentally observed CP 
violation, there are strong hints of additional sources of CP violation beyond the phase 
in the Cabibbo-Kobayashi-Maskawa matrix (5ckm)- The strongest motivation for this 
suggestion is that the strength of CP violation in SM is not sufficient to explain the 
cosmological baryon asymmetry of our universe 

In supersymmetric (SUSY) extensions of the SM, there are new CP violating phases 
which arise from the complexity of the soft SUSY breaking terms and the SUSY preserving 
//-parameter. These new phases have significant implications and can modify the SM 
predictions in CP violating phenomena. In particular, they would give large contributions 
to the electric dipole moment (EDM) of the neutron and mercury atom ||, to CP violating 
parameters (e and e') of K — K system and to the CP asymmetries in the B — B 

system These phases can be classified into two categories. The first category 

includes flavor-independent phases such as the phases of the //-parameter, £>-parameter, 
gaugino masses and the overall phase of the trilinear couplings. The other category 
includes the flavor- dependent phases, i.e. the phases of the off-diagonal elements of Aij 
and phases in the squark mass matrix raf,-. Two of the flavor- independent phases can be 
eliminated by the U(1)r and U(1)pq transformations. 

However, the non-observation of EDMs imposes a stringent constraint on flavor- 
independent SUSY phases, the so-called SUSY CP problem. Four solutions to this 
problem have been indicated so far. In the first, the CP is an approximate symmetry 
and therefore the CP violating phases are very small (~ 10~ 2 ). However, the large CP 
asymmetry in B — > ipK s decay observed by BaBar |T2J and Belle |T3j implies that CP 
is significantly violated and disfavors this possibility. In the second, the sfermion of the 
first two generations are very heavy while the third generation remains light. However, 
it was shown that in order to satisfy the EDM of the mercury atom the sfermion 
masses have to be of order 10 TeV, which leads to a large hierarchy between SUSY and 
electroweak scales. A third possibility is that the full computation of SUSY contributions 
to the EDMs involves accidental cancellations among various contributions which may 
allow for regions of parameter space with large phases. However, it has been found 
that such EDM cancellation can not occur simultaneously for the electron, neutron and 
mercury. Finally, the SUSY CP violation can have a flavor off-diagonal character just as 
in the SM |b|. In this latter possibility, the origin of the CP violation is assumed to be 
closely related to the origin of flavor structures rather than the origin of SUSY breaking. 
Thus the flavor blind quantities as the //-term and gaugino masses are real. 
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The class of models with flavor-dependent CP violation are based on two major as- 
sumptions. First, the SUSY breaking does not break CP, which can happen in some 



explicit string models Second, the flavor structure of the model is Hermitian, i.e., 
the Yukawa matrices and the A-terms have to be Hermitian, in order to ensure that the 
diagonal elements of the A-terms are real in any basis and do not induce unacceptably 
large EDMs. 

In Ref. [ 14j , it was shown that within these assumptions, the supersymmetric CP 



problem can naturally be resolved and a correlation between the CP asymmetry of the 
B — > Xg'-f decay and the EDMs is predicted. However, it was also found that SUSY 
contribution to Ek is, in general, very small and also that the dominant gluino contribution 
to e'/e is negligible due to a cancellation between the contributions involving (Sf 2 )LR and 
(^12) RL mass insertions. In fact, the SM prediction for Ek can be fitted with the current 
experimental data, however, due to the large uncertainties in the theoretical estimate of 
e' / e it is rather unclear if the SM result is consistent with the observed measurements by 



KTeV |L6| and NA31 [17]. Thus, it is necessary systematically to analyse the different 
SUSY contributions to e' / e to show if it is possible to saturate the observed value in SUSY 
models with Hermitian Yukawa couplings and A-terms. 

In this paper, we study more explicitly CP violation in the K and B system due to 
flavor dependent phases in SUSY models with the Hermiticity assumption. We show that, 
in order to saturate Ek in a viable region of parameter space non-universality between 
the squark soft masses is required. This non-universality is also essential to enhance the 
chargino contribution to e'/e. We demonstrate that the SUSY flavor off-diagonal phases 
have significant implications on the direct CP asymmetry in the B — > X57 decay, while 
their effect on the CP asymmetry in the B — ► ipKs decay are negligible. 

The paper is organized in the following way. In section 2 we start by emphasizing the 
possibility of obtaining Hermitian Yukawa couplings in SUSY models with SU(3) flavor 
symmetry, and then show that the EDMs in this class of models are one or two order 
of magnitude below the experimental constraints without any fine-tunning. Section 3 is 
devoted to the study of CP violation in the Kaon system. In section 4 we consider the 
CP violation in the 5-sector and show that in this framework the large CP asymmetry 
in the B — > ipKs decay is given by the SM contribution while the SUSY contribution is 
very small. In contrast, the SUSY contribution to the CP asymmetry in the B — > X57 
decay can be as large as ±10%. Finally, the conclusions are presented in section 5. 
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2 Hermiticity and EDM suppression 



As discussed in the introduction, an elegant solution for suppressing the EDMs in SUSY 
models is to have Hermitian flavor structures, i.e. Y a = Y at , and A a = A at . It is 
known that Hermitian Yukawa matrices can be implemented in models with left-right 
symmetry [18| and horizontal flavor symmetry In Ref.jjM], it was assumed that 



Hermitian Yukawa appeared due to a horizontal symmetry U(3)h which gets broken 
spontaneously by the VEVs of the real adjoint fields, T a , hence the Yukawa couplings are 
given as = QJL (T a ) (\ a )ij, where gn is a coupling constant of order one, M is a mass 
scale much higher than the electroweak scale, A a , for a = 1, ..,8 are Gell-Mann matrices 
and A is proportional to the unit matrix. However, the real fields T a may only arise from 
non-supersymmetric sector. In fact, if T a are the scalar components of chiral multiplets, 
they are intrinsically complex and Hermitian Yukawas arise only if the VEVs of T a are 
real. 

As we will show in the following subsection, it is possible to obtain Hermitian Yukawa 
couplings in SUSY models with SU (3) flavor symmetry broken by complex VEVs of scalar 
fields 4> a and 0° in the triplet and antitriplet representation of SU(3) respectively. The 
local SU(3) flavor symmetry provides a dynamical origin for the observed fermion masses 
and a natural explanation for three quark-lepton families. A considerable amount of work 
has been done concerning the implications of this symmetry on solving the fermion flavor 



problem [20 



2.1 Hermitian Yukawa from SU(3) flavor symmetry 

Here, we show that Hermitian Yukawa couplings can be motivated by super symmetric 
models with flavor symmetry SU(3)p- We consider a SUSY model with the gauge group 
Gsm x SU(3)p, where Gsm refers to the standard model gauge group. Under SU(3)p, 
the matter content of the MSSM is assigned the following quantum numbers: 

{Q a ,L a } = 3 and {<, d c a , e c a } = 3, (1) 

while the MSSM Higgs are singlets under the SU(3)f and have the charges {H u , H d } = 1. 
The extra Higgs fields that are used to break SU(3)f are (p a = 3 and (fi a = 3, a = 1, 2, 3. 

In this class of models, the lowest dimensional SU(3)p invariant operators in the 
superpotential, which are responsible for generating the fermion masses, are given by 

W^Yuk = h u Q a u c b H u ^ + h d Q a d c b H d ^ + h { L a e c h E d ^-. (2) 



Thus for \^rj = {(p) , i.e., (<j) a ) = v a e lipa and i^pbj = v b e t(fb , the Yukawa couplings are 
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(3) 



given by 

Y u - h VaVb A>P«-<Pb) 

1 ab n u M2 e 

Similar expressions hold for Y d and Y l . Eq.(|3]) clearly displays the usual form for Hermi- 
tian Yukawa couplings. Now let us discuss briefly the minimization of the scalar potential 
of the fields. The most general renormalizable superpotential involving these Higgs 
fields has the form 

w = v<p a r + \<p a rs + w(s) , (4) 

where the S field is a singlet under both Gsm and SU(3)p. The requirement that the 0's 
and S fields do not contribute to SUSY breaking implies that F<p = Fs = 0. The scalar 
potential is given by 
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V 



/i0 a + A0 a S + |/i0 a + \<j> a S\' 



A0 a a + 



dW 



dS 



+ V D + V SB 



where 



/ 6 



(5) 



(6) 



Note that <ft a is an anti-triplet under SU(3) p as well as a so the above potential is SU (3) p 
invariant. In the above equation, T b are the generators of the SU(3)p group, and the sum 
extends over all these generators. Finally we assume the following soft SUSY breaking 
terms 



Vf 



SB 



m 



0a I 



m- 



+ 



^ a S + B d 



+ h.c. 



(7) 



The minimization of the scalar potential with respect to <p a and a depends on the soft 
SUSY breaking terms and, for a particular choice of these parameters, one can obtain the 
following VEV's for a and 0& 



*a) = v a e 



'■fa 



-np b 



(8) 



as required in order to get Hermitian Yukawa textures. 

Furthermore, since 4> a 4> a is a singlet under both the Gsm and SU(3)p, it can couple to 
H u and Hd to generate the /x-term. In this case we can have the following leading term 
in the superpotential: 



M 



-H„Hd- 



(9) 



Thus, the //-term will be given, after the SU(3)p is spontaneously broken, by /t = v 2 a jM 
which is real. However the SU(3)p symmetry, like any other flavor symmetry and left- 
right symmetry, can not guarantee the reality of the gaugino masses and we have to make 
the additional assumption that the SUSY breaking dynamics conserves CP. This seems 



natural if CP breaking is associated with the origin of the flavor structure |H| 
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2.2 EDM-free SUSY CP violating phases 



We have shown that Hermitian Yukawa matrices and a real //-term can arise naturally in 
models with SU(3) p symmetry and the assumption that CP violation and SUSY breaking 
have different origins leads to arg(Mj) = where Mj are the gaugino masses. In this case 



the 1-terms can also be Hermitian and the EDM problem is naturally avoided 14, 18 



In supergravity models, the trilinear parameters are given in terms of the Kahler 
potential and the Yukawa couplings 



A afh = F m [K m + d m In Y a/3l - d m \n{K a K K n 



(10) 



where the Latin indices refer to the hidden sector fields that break SUSY and the Greek 
indices refer to the observable fields. According to our previous assumption, the F m is 
real. Also K a and K m are always real, thus the 1-terms are Hermitian if the derivatives of 
the Kahler potential are either generation-independent or the same for the left and right 
fields of the same generation, i.e., if Kq l K Ur = Kq l K Ur . These conditions are usually 
satisfied in string models. It is interesting to note that although the SUSY breaking does 
not bring in new source of CP violating, the trilinear soft parameters involve off-diagonal 
CP violating phases of 0(1). This stems from the contribution of the term <9 m In Yq,^, 



which has been found to be significant and sometimes even dominant in string models [21 



In what follows, we will show that these phases are unconstrained by the EDMs and will 
study their phenomenological implications in the K and B systems. 

The relevant quantities appearing in the soft Lagrangian are (Y* 4 )^- = (Y g ) i j(A q ) i j 
(indices not summed) which are also Hermitian at the GUT scale. For the sake of defi- 
niteness, we consider the following Hermitian Yukawa matrices at the GUT scale 



/ 5.94 x 10" 



Y" 



-10" 



10~ 3 i 
5.07 x 10- 3 



V -2.03 x 10- 2 -2.03 x 10- 



V 



6.84 x 10~ 3 
1.05 - 0.947 i) x 10~ 2 
-0.023 



1.05 



-2.03 x 10~ 2 \ 
2.03 x 10~ 5 i 

1 J 
- 0.947 i) x 10~ 2 

0.0489 
-0.0368 i 



-0.023 \ 
0.0368 % 

1 / 



(11) 



(12) 



where \ u = m t /v sin/3 and = rrib/v cos (3. These matrices reproduce, at low energy, the 
quark masses and the CKM matrix. The renormalization group (RG) evolution of Yukawa 
couplings and the A terms slightly violate the Hermiticity. Therefore, the resulting Y q A at 
the electroweak scale has very small non-zero phases in the diagonal elements (due to the 
large suppression from the off-diagonal entries of the Yukawa). However, what matters is 
the relevant phases appearing in the squark mass insertions in the super-CKM basis, i.e., 
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the basis where the Yukawa matrices are diagonalized by a unitary transformation of the 
quark superfields Ul,r and Dl,r (Note that since the Yukawas are Hermitian matrices 
they are diagonalized by one unitary transformation,^. e.Vl = V#): 

Ul,r -> V u U L , R7 Dl,r -> V d Dl,r-, 

yu yuT yu yu* = diag ^ ^ ^ 

Y d _> V dT Y d V d * = diag(h d ,h s ,h b ). (13) 

Accordingly the trilinear terms transform as Y q A — > V qT Y q A V q * . Thus the Y q A stay 
Hermitian to a very good degree in the super-CKM basis. Therefore, the imaginary parts 
of the flavor conserving mass insertions 

{S$ u) )lr = — 2 [(V^YfV*) .. v m - ^% 2 (1)} , (14) 
Trig L x %% 

that appear in the EDM calculations are suppressed. In the above formula the rriq refers 
to the average squark mass and Vi = (H?) /y/2. 

The effective Hamiltonian for the EDM of a fermion / containing dimension-5 and 6 
operators is given by |22| 



Hfi )M = J2C l (fi)O l + h.c, (15) 



where Oi are given by 



Oi = -\fo^fF^ 2 = - l -fa^fG a ^ 3 = - l -f abc Gl p G b SG c x<T e^. (16) 

Oi t 2 refer to the electric and chromoelectric dipole moment operators and 0% to the 
Weinberg three gluino operator. All these operators can contribute to the quark EDM 
while only 0\ contributes to the electron EDM, i.e., the Wilson coefficients Cf and Cf 
of the electron are identically zero. The supersymmetric contributions to the Wilson 
coefficients of the quark result from the 1-loop gluino, chargino, and neutralino exchange 
diagrams and also the 2-loop gluino-quark-squark diagram. As emphasized in Ref. ||, 
the most stringent constraint on the SUSY CP violating phases comes from the recent 
experimental bounds on the EDMs of the neutron and mercury atom. Therefore we will 
not discuss the electron EDM here. 

The EDMs of quarks, using the naive dimension analysis, are given by 

e eA 

^ = ^lCi + - 7/2 C 2 + — 773 C 3 , (17) 
47r Ait 

where the QCD correction factors are 7/1 = 1.53, 7/ 2 — t/ 3 ~ 3.4, and A ~ 1.19 GeV is the 
chiral symmetry breaking scale. The dominant 1-loop gluino contributions to the Wilson 
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coefficients of the down and up quarks are given by 



2 ea s 

3 V 



mi 



1 g s a s m § f <-d(u)s , •> 



7T m~ 



(18) 
(19) 



Here nig is the gluino mass and the function M\(x) is defined by 

1 + 4x — 5x 2 + 4x ln(x) + 2a; 2 ln(x) 



Afi(z) 
M 2 (x) 



2(1 -x) 4 

22 - 20x - 2x 2 + 9 ln(x) + 16x ln(x) - x 2 ln(x) 
3(1 -xY 



(20) 
(21) 



with x = ni~/m^. In the quark model, the EDM of the neutron is given by d n = ^(Add—d u ) 
and the current experimental bound |23] 

d n < 6.3 x 1(T 26 e cm (90%C.L.) 



(22) 



leads to the constraint lm(5fi )lr ^ 10 6 — 10 7 . However, it turns out that the recent 



experimental limit on the EDM of the mercury atom [El 



d Ho < 2.1 x 10~ 28 e cm 



(23) 



implies stronger bounds on these mass insertions (more than an order of magnitude 
more stringent than those imposed by the EDM of the neutron) and in addition to 
110.(822) lr % 10~ 5 — 10~ 6 0, due considerable contributions from the strange quark to 
the mercury EDM. Recall that the mercury EDM is sensitive to the chromoelectric EDM 
of quarks (Cf) and the limit in Eq. (p3|) can be translated into |C| — C 2 — 0.012Cf |/g s < 
7 x 10~ 27 cm ||. 

We start our analysis by revisiting the EDM constraints on the flavor off diagonal 
phases of SUSY models with Hermitian Yukawa as in Eq.flTTD and the following Hermitian 

A-terms: 



Ad — A,. 



( A n A 12 e^ 12 A 13 \ 
An e"^ 12 A 22 A 23 e^ 23 



V A 13 e-^ A 23 e 



V-p'2-3 



A 



(24) 



/ 



We also assume that the soft scalar masses and gaugino masses M a are given by 



M a = mi/2, a = 1,2, 3, 

m Q = = m H 2 = m l, 

m u = m D = m o diag{l, 81, 6 2 }, 



(25) 
(26) 
(27) 
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where the parameters Si and Aij can vary in the ranges [0, 1] and [—3, 3] respectively. 
Note that in most string inspired models, the squark mass matrices are diagonal but not 
necessary universal. The non-universality of the squark masses is not constrained by the 
EDMs. However, this non-universality (specially between the first two generations of the 
squark doublets) is severely constrained by AMk and Ek- 

In Fig. 1 we display scatter plots for the neutron and mercury EDMs versus the 
phase y?i2 for tan/3 = 5, m = m 1 / 2 = 200 GeV, are scanned in the range [—3,3], 
and the phases <pi3 and ^23 are randomly selected in the range [0,7r]. As stated above, 
the parameters Si are irrelevant for the EDM calculations and we set them here to one. 
Finally, since \x is real the EDM results display very little dependence on tan /3. 

It is important to mention that we have also imposed the constraints which come 
from the requirement of absence of charge and colour breaking minima as well as the 



requirement that the scalar potential be bounded from below [2B|. These conditions may 
be automatically satisfied in minimal SUSY models, however in models with non-universal 
A-terms they have to be explicitly checked. In fact, sometimes these constraints are even 
stronger than the usual bounds set by the flavor changing neutral currents [E7| . 

As can be seen from Fig. 1, the EDMs do not exceed the experimental bounds for 
most of the parameter space. Generally, they are one or two order of magnitude below the 
present limit, and the flavor-off diagonal phases of the A-terms can be 0(1) without fine- 
tunning. The points that lead to mercury EDM above the experimental bound correspond 
to (f23 — tt/2. This phase induces a considerable contribution to the chromoelectric EDM 
of the strange quark Cf . Thus the compatibility with mercury EDM experiment requires 
that the phase ^23 should be slightly smaller than tt/2. 
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3 CP violation in the Kaon system 



We have shown in the previous section that the Hermiticity of the Yukawa couplings and 
A-terms allows the existence of large off-diagonal SUSY CP violating phases while keeping 
the EDMs sufficiently small. However, the important question to address is whether these 
"EDM-free" phases can have any implication on other CP violation experiments. In this 
section, we will concentrate on possible effects in the kaon system. 

Recently, it has been pointed out that, in SUSY models with generic non-degenerate 
A-terms (where the phases of the diagonal elements are set to be very small by hand in 
any basis to satisfy the EDM bounds), it is possible to have large effects in CP violation 
observables, in particular ek and e'/e [||-|5|]. However, as we will show below, in the 
Hermitian scenario the situation is quite different and it is not straightforward to realize 
the above mentioned mechanism (which relies on the LR down squark mass insertions) 
to obtain significant SUSY contribution to Ek and e'/e. As shown in Ref . []l4f , the typical 
values of ex in this class of models are smaller than the experimental measurement (even 
with very large off-diagonal elements, Aij ~ 5mo). Moreover, it turns out that e'/e is also 
very small (~ lCT 6 ) due to a severe cancellation between the different contributions. In 
the following, we will show that we have to consider the other SUSY contributions from 
the LL and RR sectors in order to saturate both ek and e'/e. 

3.1 Indirect CP violation 

In the kaon system and due to a CP violation in K° — K° mixing, the neutral kaon mass 
eigenstates are superpositions of CP-even (Kg) and CP-odd (Kl) components. However, 
the CP-odd Kl decays into two pions through its small CP-even component. This decay, 
Kl — > 7T7T, was the the first observation of CP violation. The measure for the indirect CP 
violation is defined as 

A(K L -> 7T7T) 

£K ~ A(K S - 7T7T) • (28) 

The experimental value for this parameter is Ek — 2.28 x 10~ 3 . Generally, Ek can be 
calculated via 

E * = v^ lm{K ° lH ^ v< ° } ■ (29) 

Here H^f =2 is the effective Hamiltonian for the AS = 2 transition. It can be expressed 
via the Operator Product Expansion as 

Hg = * = Y i C i (ji)Q i + h.c. , (30) 

i 

where Ci(p) are the Wilson coefficients and Qi are the relevant local operators, which 
are given in Ref. |[29|| . The main uncertainty in this calculation arises from the matrix 
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elements of Qi, whereas the Wilson coefficients can be reliably calculated at high energies 
and evolved down to low energies via the RG running. 

The K° — K° transition can be generated through the box diagrams with W, Higgs, 
neutralino, gluino, and chargino exchange. The off-diagonal entry in the kaon mass 
matrix, M u = (K°\H^ S=2 \K°), is given by 

M 12 = M^ + M^ 2 + Mf 2 +M\ 2 +M£. (31) 



The SM contribution can be written as 130 



GlM, 



2 



m ,m = r±_fL F ^M K B K F* . (32) 

For the specific Hermitian Yukawa ansatz we are considering, we find the following SM 
contributions: 

e| M ~ 1.8 x 1(T 3 . (33) 

We see that the SM prediction for ek is close to the measured value. However, a precise 
prediction cannot be made due to the hadronic and CKM uncertainties. 

Now let us turn to the supersymmetric contributions. The Higgs and the neutralino 
contributions are very suppressed and can be neglected. The chargino contribution to the 
K° — K° mixing is given by 

M & = TtJl \M K f K BM (j2K* a2 (8t L ) ab K bl ) £ iVam^Hix^xj), (34) 

where X{ = (m-+/mg) 2 , K refers the CKM matrix, a, b are the flavor indices, i,j label the 
chargino mass eigenstates, and V is the matrix that is used for diagonalizing the chargino 



mass matrix. The loop function H(xi,Xj) is given in Ref. |pl||. However, the chargino 
contribution can be significant only if there is a large LL mixing in the up- sector, namely 
Im(5^ L ) 2 i ~ 1CT 3 and Re(5^ L ) 2 i ~ 1CT 2 [^TJ. Such mixing can not be accommodated 
with the universal scalar masses assumption (i.e., 5i = 1). In this case, the values of 
the Im((5^ i ) are of order 10~ 6 which leads to a negligible chargino contribution to Ek- 
With non-universal soft scalar masses (5i ^ 1) a possible enhancement in the chargino 
contribution is expected however, as we will show, this non-universality also leads to a 
larger enhancement in the gluino contribution. So, the dominant SUSY contribution to 
the K—K mixing in this class of models will be provided by the gluino exchange diagrams. 
The gluino contribution to the AS* = 2 effective Hamiltonian is given by ]29| 



H ^ S=2 = ^^i^M 2A ^ x ^ + QQ QM X )) + (^ 2 )L(24gi*/ 6 (*) 

ii 



216m? 



+66g 1 / 6 (a;)) + (^ 2 ) iL (^ 2 ) M (504Q 4 ^/ 6 (x) - 72Q 4 f 6 (x) + 2AQ 5 xf 6 (x) 
+120Q 5 f 6 {x)) + (5i 2 f RL (204Q 2 xf 6 {x) - 36Q 3 x/ 6 (x)) + (6 d 12 ) 2 LR (204Q 2 xf 6 (x) 
-36Q 3 x/ 6 (x)) + (S d 2 ) LR (5 d 12 ) RL (-132Q 4 / 6 (x) - 180Q 5 / 6 (x)) }, (35) 

where x = m~/m|, is the average squark mass, is the gluino mass, and the functions 
fa(x), h(x) are given by 

, , s 6(1 + 3x) In x + x 3 - 9x 2 - 9x + 17 

/■(*) = ' ( 36 ) 



/eW = — TVs • ( 3T ) 



6(1 + 3x) In x + x 3 - 9x 2 - 9x + 17 
6(x- 1) 

The matrix elements of the operators Qi between the .fT-meson states in the vacuum 
insertion approximation (VIA), where B — 1, can be found in Ref . ||28|| . The VIA generally 
gives only a rough estimate, so other methods, e.g. lattice QCD, are required to obtain a 
more realistic value. The matrix elements of the renormalized operators can be written 

as 



29 



(j^lQiMI* ) = ^k&B^I (38) 



2 



(K°\QM\K°) = -A ( M « j MKf 2 B2{ l m 
(K°\QM\K°) - -sj ^J 'w (40) 

<m</0IJO = -a UM+'^J '^^w- (41) 

where Qi{jj) are the operators renormalized at the scale /i. The expressions of the matrix 



elements of the operators Qi_3 are valid for the operators Qi_3 P5fl , and for /x = 2 GeV 
we have ||29fl 



= 0.60, B 2 (^) = 0.66, B 3 (fi) = 1.05, 5 4 (/i) = 1.03, B b (jj.) = 0.73. (42) 

Using these values, the gluino contribution to the K° — K° can be calculated via Eq . (|29|) . 
As mentioned above, for universal soft scalar masses the LL and RR insertions are gener- 
ated only through the RG running and can be neglected. The LR and RL mass insertions 
appear at the tree level and may have tangible effects. It is worth mentioning that, the 
RL and LR mass insertions contribute with the same sign in Eq.(^) and for Hermitian 
A-terms they are almost equal, so no cancellation between these two contributions occurs. 

In Fig. 2 we plot the values of \br\ versus the phase tpx 2 for Si = 1 and the other 
parameters are chosen as in Fig. 1. From this figure, we conclude that the SUSY contri- 
bution with Hermitian Yukawa and universal soft scalar masses can not account for the 
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Figure 2: The gluino contribution to \ek\ as a function of <p\2 for Si = 1, tan/3 = 5, and 
mo = mi/2 = 200 GeV. 



experimentally observed indirect CP violation in the Kaon system. In Ref . ||14j| values for 



Sk ~ 10~ 3 have been obtained but these values require light gaugino mass (mi/ 2 ~ 100 
GeV) which is now excluded by the new experimental limits on the mass of the lightets 
Higgs. Also it requires that the magnitude of the off-diagonal entries of the v4-terms 
should be much larger (at least five times larger) than the diagonal ones, which looks 
unnatural. 

A possible way to enhance Ek is to have non-universal soft squark masses at GUT 
scale. As mentioned above, the soft scalar masses are not necessarily universal in generic 
SUSY models and their non-universality is not constrained by the EDMs. Thus for Si ^ 1 
the mass insertion (5f 2 )RR is enhanced and we can easily saturate Ek through the gluino 
contribution. To see this more explicitly, let us consider the LL and RR squark mass 
matrices in the super-CKM basis 

(Mj) iL ~ U d t M 2 Q V d , 

(Mf) RR ~ V d HM 2 D ) T V d . (43) 

Due to the universality assumption of Mq at GUT scale, the matrix [Mfj remains ap- 
proximately universal and the mass insertions (Sf 2 )LL are sufficiently small (lm(5f 2 )LL ~ 
10~ 5 ). However, since the masses of the squark singlets M|> are not universal, Eq . (|27|) . 
sizeable off-diagonal elements in (Mj) are obtained. We find that Re(5f 2 )i?R is en- 
hanced from ~ 10 -7 in the universal case (Si = 1) to ~ 10~ 3 for Si ~ 0.7 while the 
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Figure 3: The value of \ek\ as a function of the parameters 5% (upper curve) and 5 2 (lower 
curve) for tan/3 = 5, and mo = "^1/2 = 200 GeV. 



imagenary part remains the same, of order 10 . Thus, in this case, we have 



Im ((5f 2 ) LL (Si 2 ) RR ) I ~ ^/\Re(5 d l2 ) RR lm(5f 2 ) LL \ ~ 10" 



which is the required value in order to saturate the observed result of Ek |28| . 

In Fig. 3 we show the dependence of \ek\ on the parameters Si. There, the two curves, 
from top to bottom, correspond to the values of \ek\ versus 8\ (for 5 2 = 1) and S 2 (for 
5i — 1) respectively. As explained above, in this scenario the main contribution of ek is 
due to LL and RR sectors and the LR sector has essentially no effect on e k . We also see 
that any non-universality between the soft scalar masses of the third and the first two 
generations can not lead to significant contribution to Ek and some spliting between the 
scalar masses of the first two generations is necessary. This stems from the fact that the 
effect of the third generation on the mass insertion (Sf 2 ) RR is suppressed by V13 ~ C(10~ 2 ) 
while V12 ~ sin 6c- Finally, as we can see from this figure, in order to avoid over saturation 
of the experimental value of Ek, the parameter S\ should be of order 0.8. 

3.2 Direct CP violation 

Next let us consider SUSY contributions to e' / e. The ratio e' / e is a measure of direct CP 
violation in the K — > tttt decays and is given by 

£ 'l £ = ~ /O M 5 A ( lmA ° ~ ~ lmA >) ' ^ 

V2 \e\ ReA \ uj J 

where A^^ 2 are the amplitudes for the AI = 1/2,3/2 transitions, and u = ReA 2 /ReAo ~ 
1/22 reflects the strong enhancement of AI = 1/2 transitions over those with AI = 3/2. 
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Experimentally it has been found to be Ke(e'/e) ~ 1.9 x 10~ 3 which provides firm evidence 
for the existence of direct CP violation. The Imv4 0i 2 are calculated from the general low 
energy effective Hamiltonian for AS = 1 transition, 

H* s=1 = J2Q(riO l + h.c., (45) 

i 

where Cj are the Wilson coefficients and the list of the relevant operators Oj for this 
transition is given in Ref. |32- |34fl . Let us recall here that these operators can be classified 



into three categories. The first category includes dimension six operators: which 
refer to the current-current operators, for QCD penguin operators and O7-10 for 

electroweak penguin operators ||32|| . The second category includes dimension five opera- 
tors: magnetic- and electric-dipole penguin operators O g and 7 which are induced by the 



gluino exchange | 33fl . The third category includes the only dimension four operator Oz 



generated by the sdZ vertex which is mediated by chargino exchanges [[34]. In addition, 
one should take into account Oi operators which are obtained from Oj by the exchange 
L <-> R. 

In spite of the presence of this large number of operators that in principle can con- 
tribute to s'/e, it is remarkable that few of them can give significant contributions. As 
we will discuss below, this is due to the suppression of the matrix elements and/or the 
associated Wilson coefficients of most of the operators. The SM contribution to e'/e is 
dominated by the operators Q G and Q$, and can be expressed as 

Re (^lEllMU, (46) 

where Aj = V* s Vid and the function F E i is given in Ref.@]. By using our Hermitian 
Yukawa in Eq.([ll|) we get 

e'/e ~ 7.5 x 10" 4 . (47) 

Again, the SM prediction is below the observed value. Nevertheless, this estimat can 
not be considered as a firm conclusion for a new physics beyond the SM since there are 
significant hadronic uncertainties are involved. 

The supersymmetric contribution to e' / e depends on the flavor structure of the SUSY 
model. It is known that, in a minimal flavor SUSY model, it is not possible to generate a 
sizeable contribution to e'/e even if the SUSY phases are assumed to be large. Recently, 
it has been pointed out that with non-degenerate A-terms the gluino contribution to e'/ e 
can naturally be enhanced to saturate the observed value [f§,0. Indeed, in this scenario, 
the LR mass insertions can have large imaginary parts and the chromomagnetic operator 
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Og gives the dominant contribution to e'/e 



Re 



11 V3 



m. 



n 



647r 2 |e|Re74 m s + F% 



m K m T 



Im 



(48) 



where C g is the Wilson coefficient associated with the operator O g , given by 



Co 



m\ 



1 

3' 



(Sf 2 )LL [ ~M 3 (x) - 3M 4 (x] 



TTlg 



--M 1 {x)-3M 2 {x] 



(49) 



where the functions Mi(x) can be found in Ref . 1 28fl and x = my mi 



Using these relations, one finds that in order to saturate e'/e from the gluino contribu- 
tion the imaginary parts of the LR mass insertions for x ~ 1 should satisfy \m(5f 2 ) lr ~ 
10~ 5 . Such values can easily be obtained in this class of models. However, as mentioned 
above, in the case of Hermitian A-terms and Yukawa couplings we have (Sf 2 ) lr ~ (Sf 2 )RL, 



hence we get Im 



C n Cr, 



Im 



hi) RL 



10 6 which leads to a negligible 



gluino contribution to e'/e [Q. It is worth noticing that, due to the universality assump- 
tion of M$ l , the imaginary part of the mass insertion (5f 2 )z,L is of order 10 5 . So its 
contribution to C g is negligible with respect to the LR one which is enhanced by the ratio 
mg/m s . To achieve the required contribution to e'/e from the LL sector, one has to relax 
this universality assumption to get Im(5f 2 )iL ~ 10~ 2 . However, as we will discuss below, 
in this case the chargino contribution is also enhanced and becomes dominant. 

Now we turn to the chargino contributions. The dominant contribution is found to be 



due to the terms proportional to a single mass insertion 

= lm\Y^K2(K b )LLK bl \ /;•!,, 




QXJ 



(50) 



a, b 



The function F E /(x qx ), where x qx 



m\±/m\ 



is given in [plj] . We find that the contri- 
butions involving a double mass insertion, like those arising from the supersymmetric 
effective sdZ, can not give any significant contribution, however we take them into ac- 
count. The above contribution is dominated by (5" 2 )ll an d in order to account for e'/e 
entirely from the chargino exchange the up sector has to employ a large LL mixing. Again, 
with universal Mq, Im(5" 2 ).LL ~ 10~ 6 and the chargino contributions (as the gluino one) 
to e'/e is negligible. 

Finally, we consider another possibility proposed by Kagan and Neubert to obtain a 
large contribution to e'/ e [[35] . It is important to note that in the previous mechanisms to 
generate e'/e one is tacitly assuming that AI =1/2 transitions are dominant and that the 
AI = 3/2 ones are suppressed as in the SM. However, in Ref. |35| it was shown that it is 
possible to generate a large e'/e from the AI = 3/2 penguin operators. This mechanism 
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relies on the LL mass insertion {S 21 )ll and requires isospin violation in the squark masses 
(rna 7^ mj). In this case, the relevant AS* = 1 gluino box diagrams lead to the effective 



Hamiltonian 135 



#cfr = %E (Ci(ji)Qi + CA) + h.c. (51) 

V 2 i=3 



where 




Ql — (d a S a ) V -A {qpQp)v+A, Q2 — (d a S(3)v-A (Qf3S a ) V +A, (52) 

Q 3 = (d a s a ) v ^ A (g /3 g ( g)y_ J 4, Q 4 = (d a sp) V -A (qpq a )v-A, (53) 

and the operators Qi are obtained from Qi by exchanging L <-> i?. It turns out that the 
SUSY A/ = 3/2 contribution to e'/e is given by |35| 

l^f^^Vs^C^^^M^W (54) 

V m 9 J 

Here, = x L d ' R = (^f^\ Bf\m c ) ~ 1 and the function K(x,y,z) is 

given by 

32 2 

^(^, y,z)^ — [f(x, y) - /(a;, + — [g(x, y) - g(x, z)] , (55) 

where f(x,y) and g(x,y) are given in Ref . ||35|| . It is clear from the above equation that 
for m^ R = rriu R , the function K(x^,x^,x^) vanishes identically. Thus, a mass splitting 
between the right-handed squark mass is necessary to get large contributions to e'/e 
through this mechanism. Furthermore, the lm(Sf 2 )LL has to be of order O(10~ 3 — 10 -2 ) 
to saturate the observed value of e'/e. It is clear that, with universal Mq, this contribution 
can not give any significant value for e'/e. 

Now we relax the universality assumption of M? at GUT scale to enhance the mass 
insertion {Sf 2 )LL and saturate the experimental value of e'/e. As mentioned above, the 
non-universality between the first two generation of M~ is very constrained by AMk and 
ex- Therefore we just assume that the mass of third generation is given by S^tuq while the 
masses of the first two generations remain universal and equal to m . This deviation from 
universality provides enhancement to both ex and e'/e. We have chosen the parameters 
5i so that the total contributions of ex from chargino and gluino are consistent with the 
experimental limits. 

In Fig. 4 we present the different gluino and chargino contributions to the e'/e and 
also the total contribution versus the imaginary part of the mass insertion {Sf 2 )LL- As 
explained above, there are two sources for the gluino contributions to e'/e: the usual 
AI = 1/2 chromo magnetic dipole operator and the new AI = 3/2 penguin operators. 
Additionally, there are two sources for the chargino contribution to e'/e: the usual gluon 
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5 




and electroweak penguin diagrams with single mass insertion and the contribution from 
the SUSY effective sdZ vertex. As can be seen from this figure, the dominant contribution 
to e'/e is due to the chargino exchange with one mass insertion. It turns out that the 
chargino contribution with two mass insertions is negligible. As expected the gluino con- 
tribution via the chromomagnetic operator is also negligible due to the severe cancellation 
between the LR and RL contributions. The contribution from the AI = 3/2 operators 
does not lead to significant results for e'/e. It even becomes negative (opposite to the 
other contributions) for lm(5f 2 )LL > 2.5 x 1CT 3 . 

From this figure we conclude that a non-universality among the soft scalar masses is 
necessary to get large values of e'/e and Sk- 



4 CP violation in the B— system 

Recent results from the 5-factories have confirmed, for the first time, the existence of CP 
violation in the B-meason decays. In particular, the measurements of the CP asymmetry 



in the Bd — > ipK s decay [[121 113]] have verified that CP is significantly violated in the 



£>-sector. The time dependent CP-asymmetry (tyK s (t) * s gi yen by 
, s T(B° d (t) -> ipK s ) - T(B° d (t) -> i)K s ) 

where AMg d is the mass difference between the two mass eigenstates of the B® — B® 



system, given by AM^ d = 0.484 ±0.010 (ps) 1 f30fl . The measurements of this asymmetry 



18 



is given by 

a ^ Ks = 0.59 ±0.14 ±0.05 (BaBar) , 

a^ Ks = 0.99 ±0.14 ±0.06 (Belle). (57) 

This large CP asymmetry implies that CP is not an approximate symmetry in nature and 
that the CKM mechanism is the dominant source of CP violation. 

In the framework of the SM, the unitarity of the CKM matrix implies the following 
relation 

V ud V: b + V cd V* b + V td V t l = 0, (58) 
which can be represented as a unitarity triangle in the complex plane (p, fj) |L1| , where 

1 



and 



p = -(1 + R* - A*) , f}=y/I$-p, (59) 

b ~ iv v*i ' * ~ iv v*r ^ J 

I V crf V c6 1 l^cd^l 

The angle /3 of this unitarity triangle is related to the complex phase of element V td and 



is defined as 11 



*W= ,™-to . (61) 
' 1 p) 2 + rj 2 



The asymmetry a^K s m the SM is given by Q^k s = sin 2(3. Using the Hermitian Yukawa 
ansatz of Eq. (|TTD, we find that the SM contribution to sin 2/? is given by 

sin 2/? SM ~ 0.59. (62) 

This result is in good agreement with the most recent world average a^K a = 0.79 ± 0.10. 
Hence, any new contribution to the CP asymmetry a^K s should be very limited. In 
supersymmetric theories the AB = 2 transition, the off-diagonal element of the B d mass 
matrix, can be written as 

/ rOI ttAB=2| r0\ 

M 12 (B d ) = {BdlH cS lBd) = M™(B d ) + M^ SY (B d ). (63) 

Zm B d 

The effect of supersymmetry can be described by a dimensionless parameter r\ and a 
phase 20 d : 

2W d = M 12 (B d ) 

M?™(B d y lb4j 

where Am^ = 2\Mf^ l (B d )\r^. Thus, in the presence of SUSY contributions, the CP 
asymmetry B d — > ipK s is modified, and now we have 

a^ Ks = sin(2/3 + 26 d ) . (65) 
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Therefore, the measurement of a^K s would not determine sin 2(3 but rather sin(2/5 + 29d), 
where 

2B < = " g t 1 + Wif) ■ (66) 

From Eqs. fl5?D and ( |6"2]) we find that the allowed range for sin 28d is as follow 

- 0.22 < sin 29 d < 0.8. (67) 

In the following, we will analyse the impact of the SUSY flavor off-diagonal phases on 
the allowed values of a^K s and the implication of the above sin 29 d constraints for the 
SUSY model we are considering. The effective Hamiltonian for AB = 2 processes can be 
expressed as 

^ =2 = E«i + ^-. (68) 



where the Wilson coefficients and the relevant operators Qi can be found in Ref . p|, |T0 



The Mf 2 VSY (Bd) receives significant contributions from the gluino and chargino exchange 
box diagrams (we also included the charged Higgs contribution in our numerical analysis). 
The gluino contribution Mf 2 {Bd) can be obtained from Eq.fl35|) by changing 8f 2 8f 3 , 
M K — > M Bd , f K — > f Bd , m s — > m b . The chargino contribution is given by |36 



oi 2 , 6 2 1 

MUB d ) = iffl d B Bd ri- %l2 'M Bd ^ E ^ jki A ihj G\x^,x^,x mj ), (69) 

h,k=l i,j=l m Xj 

where the function Aijk and G'(x,y, z) can be found in Ref.[^]. In most of the param- 



eter space, we found that the chargino contribution gives the dominant effect to the CP 
asymmetry sin 26 d, while the gluino is sub-dominant. This result can be understood by 
using the mass insertion method. The gluino amplitude receives a leading contribution 
from (5i 3 )ll, since the mass insertions (<5f 3 )_L_R and (6f 3 ) B R are much smaller. However, for 
m = mi/2 ~ 200 GeV, the |(<5i 3 )ll| ~ 10~ 3 which is two orders of magnitude below the 
required value to saturate the experimental value of a^K s |l3 , so that one has negligible 



gluino contribution to sh^^. The chargino amplitude is proportional to (^i 3 )ll which 
can be enhanced by a light stop mass. 

Our results for the total SUSY contribution to the CP asymmetry sin 29 d as a function 
of the |(<5i 3 )ll| are presented in Fig. 5. Here we have also assumed m = m\/2 = 200 
GeV and tan/5 = 5. The parameters Aij, (pij and 812 are varied in their allowed regions 
fixed by the experimental limits on the EDMs, Ek and e/e. As expected, in order to 
have significant SUSY contributions to the CP asymmetry a^K s a large mixing in the LL 
sector is required (in order to enhance the dominant chargino contribution). However, 
such mixing is not allowed in our model with Hermitian Yukawas, as Fig. 5 confirmed. 



20 



0.2 



0.1 5 h 

"D 

<x> 

■tf) 



**. *% * « • «* • • 

I •* *"„• *,*•< * * # 

» * - • :i- m *.<£ . .>• . . •* • 



0.002 0.004 0.006 0.008 0.01 0.012 0.014 
I (5 d l 3 )LL I 

Figure 5: The SUSY contribution to sin 26^ and versus the (|#l 3 )ll|- 



This, in fact, is due to the Ek constraint that severely cuts off any enhancement though 
the non-universality of the soft scalar masses. 

Also as can be seen from Fig. 5, the predicted values of sin 29 d reside within the 



allowed range defined in Eq. (|67j) . Hence, in this class of models, the SM gives the leading 
contribution to the CP asymmetry a^K s , while the supersymmetric contribution is very 
small. Therefore, there is no constraint can be set on the non-universality of this class of 
models by the recent BaBar and Belle measurements. 

We conclude this section with some remarks on the CP asymmetry in B — > X s 7. 
It is known that the SM prediction for the CP asymmetry is very small, less than 1%. 
Thus, the observation of sizeable asymmetry in this decay would be a clean signal of 
new physics. The most recent result reported by the CLEO collaboration for the CP 
asymmetry in these decays is —27% < A^p 1 < 10% [0. The SUSY predictions for 
^cp 1 are strongly dependent on the flavor structure of the soft breaking terms. In the 
universal case, as in the minimal supergravity models, the asymmetry is found to be less 
than 2% ||. However, it was shown that the non-universal A-terms can result in a large 
CP-asymmetry in the B — > X s 7 || and these effects are correlated with the EDMs. 

The enhancement of A b ^p sl is due to important contributions from gluino-mediated 
diagrams, in this scenario, in addition to the usual chargino and charged Higgs contri- 
butions. The expression for the asymmetry A b ^p sl , corrected to next-to-leading order is 
given by 



A b^s~f 
A CP 



9ICVI 2 



10 



-2z {v(z)+b{z,6)) 



Im 



Im 
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+ 2 -z b(z,5)Im[C 2 C 8 s +C 2 C*]}, (70) 

where z = m 2 c /m 2 . The functions v(z) and b(z,S) and the Wilson coefficients Cj can be 
found in Ref.0,g. In the EDM-free models we are considering, we found that the flavor 
dependent phase (^23 gives a large contribution to the CP asymmetry (since the Wilson 
coefficients are proportional to {S^lr which receives a dominant contributions from A23 
entry |3J). The effect of the other flavor dependent phases on A b ^p sl is found to be very 
small. 

In our SUSY model with Hermitian A-terms, we found that the gluino contribution to 
CP asymmetry A^p 1 can be as large as 10%, which can be accessible at the B facrories. 



5 Conclusions 

In this paper, we have analysed the CP violation in supersymmetric models with Hermi- 
tian Yukawa and trilinear couplings. We emphasized that in most of the parameter space 
of this class of models the EDM of the neutron and mercury atom are two orders of mag- 
nitude smaller than the experimental limits. Furthermore, we have shown that Hermitian 
Yukawas can naturally be implemented in SUSY models with £77(3) flavor symmetry. 

We have studied the CP violation in the Kaon system. We found that in order to 
saturate Ek a small non-universality between the squark soft masses is required. We 
investigated the effects of the EDM-free, flavor off-diagonal, phases on the direct CP 
violation observable e'/e. A large SUSY contribution to this observable is possible via 
the chargino contribution. 

Finally, we considered the a^x s an d ^cp 1 CP asymmetries in 5-meason decays. We 
verified that the SM contribution to a^K s is in agreement with the recent measurements 
by BaBar and Belle, while the SUSY contributions are very small and hence no further 
constraint is imposed on the non-universality of these models. In contrast, with A b ^p sl 
the SM contribution is negligible and the SUSY contribution can be as large as ±10% 
which can be accessible at B-factories. 
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